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Abstract 

High-resolution ac-calorimetry has been carried out on dispersions of aerosils in the liquid crystal 
octyloxycyanobiphenyl (80CB) as a function of aerosil concentration and temperature spanning 
the crystal to isotropic phases. The liquid-crystal 80CB is elastically stiffer than the previously 
well studied octylcyanobiphenyl (8C-B)+aerosil system and so, general quenched random disorder 
effects and liquid-crystal specific effects can be distinguished. A double heat capacity feature 
is observed at the isotropic to nematic phase transition with an aerosil independent overlap of 
the heat capacity wings far from the transition and having a non-monotonic variation of the 
transition temperature. A crossover between low and high aerosil density behavior is observed for 
80Ci?-|-aerosil. These features are generally consistent with those on the 8C-B+aerosil system. 
Differences between these two systems in the magnitude of the transition temperature shifts, heat 
capacity suppression, and crossover aerosil density between the two regimes of behavior indicate 
a liquid crystal specific effect. The low aerosil density regime is apparently more orientationally 
disordered than the high aerosil density regime, which is more translationally disordered. An 
interpretation of these results based on a temperature dependent disorder strength is discussed. 
Finally, a detailed thermal hysteresis study has found that crystallization of a well homogenized 
sample perturbs and increases the disorder for low aerosil density samples but does not influence 
high density samples. 

PACS numbers: 64.70.Md, 61.30.Eb, 65.40.Ba 
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I. INTRODUCTION 



The effect of quenched random disorder on phase structure and transitions is an important 
area of study that continues to attract a great deal of research. Disorder is ubiquitous (ideal 
pure transitions being the exception rather than the rule in nature) and the effect on phase 
transitions can be profound. Phase transitions are modified depending on the aspect of the 
system affected by the disorder, on the dimensionality, and on the number of components to 
the order parameter. Of particular interest is a linear coupling between quenched random 
disorder and the order parameter, which allows for a random-field theoretical approach. Also, 
the order of the transition is crucial as first-order transitions have additional considerations 
compared to continuous transitions. This is due to the presence of two-phase coexistence 
(hence interfaces between ordered and disordered regions), intrinsically finite correlation 
length at the transition, and hysteresis effects for first-order compared to continuous phase 
transitions. This has made the experimental and theoretical studies of quenched random 
disorder effects at first-order transitions challenging. 

Liquid-crystals (LC) are a particularly attractive system for the study of phase transitions 
into partially ordered phases. This makes them especially interesting for the study of the 
effects of quenched random disorder (QRD), which are typically introduced by the random 
fixed dispersion of solid surfaces. In LC+aerosil systems, the quenched random disorder is 
created by a dispersed gel of aerosil particles and is varied by changing the density of aerosils 
in the dispersion. A convenient measure of the introduced disorder is the grams of silica per 



related to the 
. The aerosils 



cm 3 of liquid-crystal, denoted the conjugate silica density p$, which is direct 
surface area of solids as well as the mean-distance between solid surfaces ^, 
used are silica spheres that can hydrogen bond together to form a fractal-like random gel. 
Studies have previously been carried out by various groups on liquid crystals in an aerogel 
medium ^. Aerogels are self-supporting structures and this places a lower limit on the 
disorder strength that can be probed. By contrast, the aerosil gel provides a weaker and 
more easily controlled perturbation, and thus opens up a physically interesting regime. 

In this work, we study the effect of quenched random disorder due to a dispersed 
thixotropic aerosil gel on the weakly first-order isotropic to nematic (I-N) phase transition. 
The calorimetric results for the nematic to smectic-^4 phase transition in 80C-B+aerosil 
samples have been previously reported and were shown to be consistent with results in 
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8CB+aerosil samples 



^J. At the I-N transition, the orientational order has a finite correla- 
tion length and is established in three dimensions, which is describable by a symmetric and 
traceless 2nd rank tensor 13 (as such, it possesses only five independent components). 
Thus, nematic order belongs, in principle, to a d = 3, n = 5 Heisenberg class. However, 
by ignoring any biaxial character and aligning the orientation axis with a principle axis 
of a local frame, this tensor can be split into a scalar order parameter (S) measuring the 
magnitude of orientational order about the orientation axis and a "headless" vector called 
the nematic director n (n = —n) describing the spatial orientation of this axis. In this 
simplified view, nematic order is described on short length scales by S and on longer length 
scales by n, which is useful in describing the elastic properties of the nematic structure. 
These measures of nematic order are related to the quadrupolar nematic order parameter 
by Qij = \S(?>hihj - 5%). 

In principle, the effect of the aerosil gel network on the orientational order of the nematic 
phase is two fold. The silica gel firstly dilutes the liquid crystal and secondly creates a pre- 
ferred local orientation JfJE]]. In addition, the first-order transition from the isotropic to the 
nematic phase necessitates the formation of interfaces between coexisting domains/phases, 
which must occur within the available void spaces. The latter effect is the classic result of 
quenched random disorder; a distribution of transition temperatures due to the nucleation 
of ordered domains within voids having some size distribution. This leads to short-range 
order (SROV, a rounding of the transition, and suppression of the first-order character of the 
transition The contribution of a random preferred local orientation effect to the total 
hamiltonian can be represented as 

Wr/ = "l>(hi.fli) a (1) 

i 

where rii is the orientation of the molecules over some small region where the orientation is 
approximately constant and hi is the random influence of the silica surface. The variance 
of this random field (h 2 ) should be proportional to the density of solids dispersed in the LC 
medium. This term is squared due to the effective inversion symmetry of the molecules in 
the nematic phase. Since the nematic order parameter is quadratic in rii due to the same 
inversion symmetry p|, Eq. $1} is also linear in the order parameter and hence constitutes 
a random field (RF) interaction. Recently, Eq. ((TJ has also been interpreted as a random- 
anisotropy (RA) interaction |3] but this seems only applicable to systems describable by a 



pure vector n = 3 order parameter. The formation of interfaces and the resulting surface 
energy penalties places restrictions on the effects of (h 2 ) depending on the elasticity of 
the nematic. Light scattering measurements have shown that the nematic phase in liquid 
crystal and aerosil dispersions breaks up into large (micron size) but finite-size domains j^l- 
In addition, more extensive optical studies focussing on the nature of the nematic director 
structure well below the I-N transition have shown that the director correlation length ^ 
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These 



decays exponentially with distance, which is a hallmark of short-range order 
features are consistent with an RF interaction for nematics with QRD. 

To date, the most thoroughly studied LC+aerosil system is the dispersion of tvpe-300 

10, Q, 



x-ray scattering 
light-scattering 



aerosil in o ctylc yanobiphenyl (8CB), denoted 8CB+aerosil. Detailed calorimetric 

14 J3, x-ray intensity fluctuation spectroscopy [r| , static and dynamic 
% Ulll, and deuterium NMR [itJ studies on the nematic to smectic-y4 
(iV-SnM) and the isotropic to nematic {I-N) phase transitions of this system have shown 
that there are clear quenched random-field characteristics as well as finite-size scaling effects 

a. 

Calorimetry measurements on 8CB+aerosil samples have been particularly useful in 
yielding detailed information on both the I-N and the iV-Sm^A phase transitions Q. The 
results for both transitions show a complex dependence of the transition temperature on the 
aerosil density. While the iV-SmA heat capacity peak remains sharp and evolves towards 
3D-XY behavior with increasing silica density, the I-N behavior is more complicated. For 
silica densities below ps ~ 0.1 g cm" 3 , two heat capacity peaks, closely spaced in temper- 
ature, were observed. At higher aerosil densities, the heat capacity peaks for both the I-N 
and the iV-Smyl transitions displayed a highly smeared and non-singular features. Deu- 
terium NMR measurements on deuterated 8CB-t-aerosil dispersions, which were carried out 
over a wide range of silica densities, showed that the magnitude of the orientational order 
S below the I-N transition temperature was essentially unchanged from bulk behavior 
The amount of liquid crystal reorientation for field-cooled samples upon rotation within the 
DNMR field is small and decreases continuously with silica density up to ps = 0.094 g cm -3 
(the units will be dropped hereafter) confirming distinct low and high ps behavior. A x-ray 
intensity fluctuation spectroscopy (XIFS) study have found evidence of aerosil gel dynamics 
in 8CB+aerosil dispersions indicating an elastic coupling between the gel and LC 16j . The 
optical, calorimetric, DNMR, and XIFS results all appear to be consistent with a model in 
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which director fluctuations are suppressed with increasing aerosil density. 

The present work focusses on a different liquid crystal - octyloxycyanobiphenyl (80CB) 
- having dispersed in it the same type of aerosil over a comparable range of silica densities as 
the well-studied 8CB+aerosil system. This liquid-crystal has several important differences 
from the closely related 8CB. The liquid-crystal 80CB has stronger smectic and nematic 
interactions than 8CB as evidenced by the higher transition temperatures, the larger bare 



correlation lengths for smectic interactions 



18| . and the larger elastic constants. More 



specifically, 80CB has a 17% larger bend, 36% larger twist, and 10% larger splay nematic 
elastic constants than 8CB (in the single elastic constant approximation, 80CB has a 
~ 20% greater than 8CB with an overall uncertainty of 5%) Thus, comparison 

of behaviors between 8C-B+aerosil and 80CB+aerosil systems allow for the isolation of 
general quenched random-disorder (QRD) effects from material specific effects, in this case 
the elasticity of the liquid-crystal host medium. 

In general, this work reveals a non-monotonic silica density dependence of the I-N and 
iV-SnM transition temperatures similar to that observed for 8C-B+aerosil but occurring over 
a larger ps range for 80CB+aerosil. The calorimetric results presented here for the I-N 
transition reveal the onset of a double transition peak for p$ > 0.1 with a ps dependence on 
the temperature distance between the two heat capacity peaks. Evidence is presented that 
the first-order character of the I-N transition continuously decreases with silica content, 
becoming approximately zero for p$ > 0.7. Over the entire range of ps studied here, the 
heat capacity temperature dependence away from the immediate vicinity of the transition 
region is bulk-like and independent of silica content. 

We speculate that the variance of the disorder (h 2 ) may change through a first-order 
transition for nematics to account for these observations. Such a variation of the disor- 
der strength may be due to the silica surfaces introducing a low-order, paranematic-like, 
boundary layer initially screening the remaining liquid crystal material. The thickness of 
this boundary layer is strongly temperature dependent in the immediate vicinity of the I-N 
transition and as it shrinks, the screening becomes weaker. 

Section [H] describes the preparation of the 80CB+aerosil dispersions as well as the ac- 
calorimetry technique employed. Given in Section IIHI is a presentation of the results. All 
results are then discussed in Section |IV| and related to results from previous LC+aerosil 
studies. Directions for future study will also be discussed. 
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II. EXPERIMENTAL TECHNIQUES 



The liquid crystal 80CB, purchased from Aldrich, was used after degassing in the 
isotropic phase for 1 hour. This liquid- crystal molecule has an aliphatic tail attached by an 
oxygen link to the rigid biphenyl core and a polar cyano head group ( M w = 307.44 g mol" 1 
). This oxygen link constitutes the sole molecular difference between 80CB and 8CB. Pure 
80CB has a weakly first-order isotropic to nematic transition at Tf N ~ 353 K and a second- 
order nematic to smectic-zl transition at T^ A ~ 340 K. At lower temperatures, the strongly 
first-order Crystal-SmA transition occurs reproducibly on heating at T^ rA ~ 328 K and, as 
usual, can be greatly supercooled. 

The hydrophilic type-300 aerosil obtained from Degussa j^j was thoroughly dried at 
~ 300°C under vacuum for a couple of hours prior to use. The hydrophilic nature of the 
aerosils arises from the hydroxyl groups covering the surface and allows the aerosil particles 
to hydrogen bond to each other. This type of bonding is weak and can be broken and 
reformed, which leads to the thixotropic nature of gels formed by aerosils in an organic 
solvent. Crystallization severely disrupts the gel. The specific surface area measured by 
the manufacturer via BET nitrogen isotherms is 300 m 2 g" 1 and each aerosil sphere is 
roughly 7-nm in diameter. However, SAXS studies have shown that the basic aerosil unit 
consists of a few of these spheres fused together during the manufacturing process Q . Each 
80C-B+aerosil sample was created by mixing appropriate quantities of liquid crystal and 
aerosil together, then dissolving the resulting mixture in spectroscopic grade (low water 
content) acetone. The resulting solution was then dispersed using an ultrasonic bath for 
about an hour. As the acetone evaporates from the mixture, a fractal-like gel forms through 
diffusion-limited aggregation. Small angle x-ray studies have shown that the aerosil gel 



dispersion has a fractal structure and no preferred orientation 111 on the micron-long length 
Jesof—o^y. 




At room temperature, 80CB is a crystalline solid even in the presence of high aerosil 
density. Care was taken to avoid crystallization of 80CB and possible damage to the 
aerosil gel, especially for low silica densities. For this calorimetry study, the mixture after 
slow solvent evaporation was allowed to crystallize, and the solid sample was transferred into 
the calorimetry cell. The cell was then sealed, the heater and thermometer attached, and 
the cell was heated into the isotropic phase. The sample was then remixed by placing the 
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assembly in an ultrasonic bath for over an hour. The cell and sample temperature was kept 
elevated during the mounting of the sealed cell into the calorimeter by maintaining current 
through the heater. This sample preparation protocol also allows a controlled entry into the 
crystal phase. However, since the cell is sealed, the in situ remix could not be inspected and 
so some small dispersion inhomogeneity may remain. 

High-resolution ac calorimetry was performed using two home-built calorimeters at WPI. 
The sample cell consisted of a silver crimped-sealed envelope ~ 10 mm long, ~ 5 mm wide, 
and ~ 0.5 mm thick (closely matching the dimensions of the heater). After the sample was 
introduced into a cell having an attached 120-0 strain-gauge heater and 1-Mfi carbon-flake 
thermistor, a constant current was placed across the heater to maintain the cell temperature 
well above Tj^. The filled cell was then placed in an ultrasonic bath to remix the sample. 
After remixing, the cell was mounted in the calorimeter, the details of which have been 



described elsewhere |2l|. In the ac-mode, power is input to the cell as P ac e lut resulting in 
temperature oscillations with amplitude T ac and a relative phase shift of (p = $ + 7r/2, where 
$ is the absolute phase shift between T ac {uj) and the input power. The specific heat at a 
heating frequency to is given by 

n _ [CfMed~ & empty] _ P ac COS(f/uJ \T ac \ — C empt y , . 

O p — — , (Z) 

^sample ^sample 

where C'^ illed and C'j iUed are the real and imaginary components of the heat capacity, C empty 
is the heat capacity of the cell and silica, m samp i e is the mass in grams of the liquid crystal 
(the total mass of the 80CB+aerosil sample was ~ 20 mg, which yielded m samp i e values in 
the range of 13 — 20 mg), and R is the thermal resistance between the cell and the bath (here, 
~ 200 K W _1 ). Equations (J2J) and © require a small correction to account for the finite 
internal thermal resistance compared to R, and this was applied to all samples studied here 



22j . Measurements were conducted at various frequencies in order to ensure the applicability 
of Eqs. (J2J) and Q by checking that C"j illed ~ through the effective iV-Sm^l transition at 
T* and that C p was independent of to. All data presented here were taken at uj = 0.1473 s _1 
at a scanning rate of less than ±100 mK h -1 , which yield essentially static C p results. All 
80C-B+aerosil samples experienced the same thermal history after mounting; six hours in 
the isotropic phase to ensure homogeneous gelation, then a slow cool deep into the smectic 



phase before beginning the first detailed scan upon heating. 



III. RESULTS 

A. General description 

The heat capacity of the pure 80CB liquid-crystal is in good agreement with previously 
published results j^J I^J. For our pure 80CB material the transition temperatures were 
Tj N = 352.53 K and T NA = 339.52 K. The I-N two-phase coexistence width was ~ 95 mK 
wide, and the A-SnM transition enthalpy was SH NA = 0.42 J g -1 . These thermal features 
indicate that the 80CB liquid-crystal used in this study was of reasonably good quality. A 
summary of the calorimetric results for pure 80CB and 80CB+aerosil samples is given in 
Table □ 

In order to determine the excess heat capacity associated with the phase transitions, 
an appropriate background was subtracted. The total sample heat capacity over a wide 
temperature range had a linear background, C p (background), subtracted to yield 

AC P — C p — C p (background) (4) 

as the excess C p due to the I-N and A-Sniv4 phase transitions. The resulting AC P data 
are shown for pure 80CB and all 80C-B+aerosil samples in Fig. over a wide temperature 
range about Tin, where the units are J K -1 per gram of liquid crystal. The transition 
temperature, Tin, is determined as the highest temperature where any nematic phase is 
present and corresponds to the highest temperature peak in Cf illed . 

As seen in Fig. ^ the AC P values away from the A-Sult4 transition and the I+N coexis- 
tence regions overlap with bulk behavior independent of silica concentration. The detailed 
variations of AC P associated with the A-Sm^ transition with p$ has been reported previ- 
ously [4]. The deviations of some of the ps = 0.05 points in the nematic phase is likely a 
consequence of sample inhomogeneity. The AC P "wings" of the I-N transition are associated 
with short-range fluctuations of nematic order. Given the simplification of the nematic order 
parameter, the short-range fluctuations in bulk nematics are mainly composed of thermal 
fluctuations of the scalar part S. For the 80CB+aerosil system, the temperature dependence 
of ACp(IN) being independent of ps suggests that thermal fluctuations of S are independent 
of disorder over the whole range of ps studied in this work. The AC P (IN) wing behavior 
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shown here for 80CB+aerosil is completelyconsistent with similar results for 8CB+aerosil 
[jj and low-density 8CB+ aerogel samples j^ . 

In stark contrast to the behavior of AC P in the one-phase regions, the two-phase coexis- 
tence region of the I-N transition exhibits strong effects of silica concentration, as shown in 
Fig. |21 From 1 K below to 0.4 K above Tin, the AC P (IN) peaks for the pure 80CB and 
the p,s = 0.036 sample are essentially the same 2(|. Upon increasing silica density, the peak 
in ACp(IN) is substantially lower in temperature relative to the peak in C'L led and consid- 
erably broader than for the bulk or the ps = 0.036 sample. See the upper panel in Fig. |21 
In addition, there is a small and very broad shoulder below the main specific heat peak, 
also seen in bulk, which moves toward the main peak with increasing p$. The nature of this 
subsidiary feature is not known. It is likely, given the similarity of the materials used here 
with the 8CB+aerosil system, that the percolation threshold for type-300 aerosil in 80CB 
is essentially the same at pp ~ 0.018 [j| and so a true gel should be present for all samples 
studied in this work. This is supported by a visual inspection of these samples holding their 
shape above the crystal melting temperature. The I-N transition regions shown in Fig. 121 
for ps ;$ 0.1 appear to be quite sensitive to small inhomogeneities in the silica gel dispersion 
and so, are grouped together. Given the dominance of a large, broad, specific heat peak and 
relatively erratic transition temperature shifts, discussed below, the effect of the silica gel 
on the I-N phase transition is strongly dependent on the quality of the dispersion. 

Beginning with the ps = 0.22 sample and for increasing silica content there is a systematic 
variation of the excess heat capacity, which is shown by the lower panel in Fig. El At 
ps = 0.22, a double heat capacity feature is observed with a sharp high-temperature peak 
ACp T corresponding closely to (but very slightly below) a sharp peak in C'j illed , followed at 
lower temperature by a broader peak ACp T also having an associated broad peak in Cf illed . 
Clearly, both are first-order signatures and they are separated by m 0.1 K. For the ps = 0.347 
and 0.489 samples, the ACp T feature remains sharp but decreases in magnitude while the 
ACp T feature becomes increasingly rounded and moves to lower temperature relative to 
ACp T by 0.15 and 0.2 K, respectively. For the ps = 0.647 sample, both heat capacity 
features are rounded and separated now by ~ 0.5 K. Over this entire range of silica density, 
the size of the C'L led peak decreased monotonically with increasing p s . Such a double I-N 
heat capacity feature was observed in 8CB+aerosil samples for silica concentrations up to 
ps ~ 0.1, but AC P exhibited a single, rounded feature above this density []]. Only a single 
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rounded AC P (IN) feature was observed for all 8 CB+ aerogel samples [2- 



B. The I-N transition enthalpies 

The I-N transition enthalpy also exhibits a dependence on aerosil concentration and 
can be a quantitative measure of the strength of the transition. For a second-order (or 
continuous) phase transition, the change in enthalpy through the transition is given by 



where the limits of integration are as wide as possible about the heat capacity peak. How- 
ever, for first-order transitions the situation is complicated by the presence of a two-phase 
coexistence region, in this work I+N, as well as a latent heat AH. The total enthalpy 
change through a first-order transition is the sum of the pretransitional enthalpy and the 
latent heat. In an ac-calorimetric measurement, AC P values observed in the two-phase re- 
gion are artificially high and frequency dependent due to partial phase conversion during a 
T ac cycle. The pretransitional enthalpy 8H is typically obtained by substituting a linearly 
truncated AC P behavior between the bounding points of the two-phase coexistence region 
into Eq. (J3J), and an independent experiment is required to determine the latent heat AH 
Q]. A direct integration of the observed AC P yields an effective transition enthalpy 5H* and 
this contains some of the latent heat contributions; thus 5H < 5H* < AH tota i = 5H + AH. 

For our analysis, the observed AC P (IN) was directly integrated over a wide temperature 
range of —25 K below to +5 K above Tin for all bulk and 80C-B+aerosil samples where the 
iV-Sim4 transition enthalpy contribution was subtracted. This will be referred to hereafter 
as the ac-enthalpy and denoted as SHj N , as it represents only a part of the total transition 
enthalpy. As seen in Fig. Q an integration of a linearly truncated AC P (IN) in the two-phase 
coexistence region over a similar range yields a pretransitional enthalpy 5Hjn = 5.13 J g _1 
that is independent of aerosil density. Integration over a similar temperature range yielded 
a pretransitional SHin value of 5.43 J g _1 for 8C-B+aerosil samples, also independent of 
silica density [jj. In addition, the integration of the imaginary heat capacity given by 
Eq. © and normalized to the LC mass, defines an imaginary transition enthalpy, referred 
to as im-enthalpy and denoted as SH'/ N , which is an indicator of the first-order character 
of the transition. Although SH'/ N is a measure of the dispersive component of the complex 




(5) 
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enthalpy, it is only approximately proportional to the transition latent heat due to the fixed-u; 
ac-technique employed in this work. As the silica content changes, the two-phase conversion 
rate may change and so alter the proportionality between 8H" N and AHj^', thus a detailed 
frequency scan for each sample would be needed to fully characterize the relationship. This 
was done for a few samples and the frequency employed in this work is sufficiently close to 
the static limit that this effect should be minimal. 

The results of both the ac- and im-enthalpy for 80CB+aerosil samples are shown in 
Fig. |3] as a function of the silica density. There is a slight variation (first increasing for 
increasing ps up to 0.220 then decreasing for larger p$) of the ac-enthalpy due mainly to 
changes in AC P values within the two-phase coexistence range since the heat capacity wings 
away from the transition are ps independent (except for p$ = 0.051, which is systematically 
high for T — Tin from —3 K to —10 K). Given the fixed-a; aspect of the technique, any 
variation observed in the ac-enthalpy in the two-phase region can be attributable to changes 
in either the dynamics or magnitude (or both) of the latent heat evolution. The small non- 
monotonic variation of SH^ N for 80CB+aerosil samples is in contrast to the systematic 
decrease of 5H* IN with increasing ps for 8CB+aerosil samples. This may reflect a difference 
in the phase conversion dynamics between 8CB and 80CB and how they are modified by 
the presence of aerosils. 

The interpretation of the im-enthalpy is more straight-forward as it is closely related to the 
latent heat of the transition. With increasing ps, the im-enthalpy appears to monotonically 
decrease to almost zero for ps = 0.647. See Fig. El This suggests that for the highest ps 
sample studied, the I-N latent heat has become nearly zero. Similar trends were observed 
for 8C-B-|-aerosil [l| where a continuous I-N transition is estimated to occur near p$ ~ 0.8. 
Also, a nearly continuous I-N transition was reported for 7CB-|-aerosil for silica densities 



near p s 



The above observations are consistent wit 



i the general view that with 



increasing QRD, first-order transitions are driven continuous 



C. Transition temperatures and crystallization 

The I-N transition temperature, defined here as the peak in C'j illed for the highest tem- 
perature feature^for the 80C-B+aerosil samples as well as those for the 8C-B+aerosil system 
taken from Ref. |l| are shown in Fig. [Has a function of silica density. For the 80CB+aerosil 
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system, Tin is essentially unchanged up to ps = 0.051 then decreases sharply by ~ 1.5 K 
at Ps — 0.078, 0.105, and 0.220. It then rises strongly for ps = 0.347, nearly recovering 
the bulk value. Upon further increase in p$, Tin decreases monotonically (with a concave 
downward character) until it is again about ~ 1.5 K below T° N for the ps = 0.647 sam- 
ple. The non-monotonic evolution of Tin with silica content for the 80CB+aerosil system 
is similar to that seen in the 8CB+aerosil system, suggesting that the initial depression 
of Tin, recovery, then continued depression is a general phenomena of quenched random 
disorder on nematics while the specific ps dependence is liquid crystal material dependent. 
Over this same range in silica density, the width of the two-phase coexistence region 5Tin 
also has a non-monotonic dependence on p$, as seen in Table However, STin is sensitive 
to local inhomogeneities of the aerosil dispersion that may account for its variation when 
ps < 0.1. Beginning at ps = 0.105, 5Tin increases monotonically by a factor of ~ 6.7 while 
Ps increases by a factor of ~ 6. The observed broadening of the two-phase coexistence width 
in nearly direct proportion with increasing QRD is generally consistent with the behavior 
of first-order transitions with quenched disorder 

The iV-Smyl pseudo-transition temperatures T* scaled by the bulk transition temperature 
T NA for 80C-B+aerosil and 8C-B+aerosil systems are shown in Fig. The pattern of 
fractional changes in the T* is essentially the same for both LC+aerosil systems with an 
initial rapid depression, recovery, then more gradual decrease with a total change of less than 
1% from T NA . The primary difference with the 8CB+aerosil system is that this behavior is 
"stretched" in ps for the 80CB+aerosil samples. This is consistent with the evolution seen 
of Tin shown in Fig. |U and described above. 

The nematic phase temperature range, ATjy = Tin~T*, normalized by the bulk nematic 
range AT N = T° N — T NA , ls shown in Fig. El While the individual transition temperature 
changes reflect the absolute stability limit of the nematic and smectic phases, ATjy reflects 
the relative stability of both phases. For 8CS+aerosil, a decrease of ~ 1% in ATjy was 
seen up to ps ~ 0.1, corresponding to the the local maximum of Tin{ps) and T*(ps). This 
was originally thought to be scatter in the data of Ref. [l]. For 80C-B+aerosil samples, 
a similar and far more pronounced 4% decrease in ATjy is seen from ps = 0.051 to 0.220. 
This decrease in ATjy reflects a greater depression of Tin than T* and indicates that in this 
range of silica density, the disorder primarily effects nematic (orientational) ordering. Upon 
further increasing ps, the nematic range begins to increase and appears to saturate at an 
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8% increase similar to that seen in the 8CB+aerosil system. This growth in the nematic 
range occurs because of the greater suppression of T* relative to Tin and so, reflects that 
above ps > 0.2, the effect of the silica gel is to mainly disorder smectic (lD-translational) 
ordering. 

The calorimetric results on the I-N transition temperature described above suggest the 
importance of sample homogeneity. As a test of the fragility of the silica gel, heat capacity 
scans were performed on a low and high density 80C-B+aerosil sample immediately be- 
fore and after crystallization of the LC. Such a thermal cycle for the ps = 0.051 sample is 
shown in Fig. and for the ps = 0.220 sample in Fig. |H1 as a function of AT = T — Tin 
in order to suppress hysteresis effects of the I-N transition. The effect of crystallization on 
the ps = 0.051 sample is striking, revealing significant distortion of the AC P signature at 
both the I-N and iV-SnM transitions. The appearance of an additional broad AC P feature 
beginning ~ 0.9 K below Tin as well as a broadened feature over the iV-Smyl transition re- 
gion after crystallization may indicate increased sample inhomogeneity, presumably caused 
by the local expulsion of silica particles as LC crystallites form. However, there are two 
puzzling aspects; (1) the shift in Tin is downward by ~ 0.7 K and (2) the specific dis- 
tortion seen in AC P about T* reveals an increased iV-Sm^l transition enthalpy. The first 
aspect is counter-intuitive as the expulsion of impurities upon crystallization should have 
moved the system closer to bulk behavior by increasing the size of pure LC domains (re- 
gions where no silica is present). The second aspect is particularly puzzling as the iV-Sm^l 
transition enthalpy after crystallization surpasses the iV-Sm^ transition enthalpy for bulk 
80CB (AH NA (before) = 0.334, AH NA (after) = 0.65, and AH NA (bulk) = 0.42 all given in 
units of J g _1 ). Visual inspection of the sample immediately after the first crystallization 
revealed no obvious inhomogeneities. A recent Raman spectroscopy study of 8CB crys- 
tallizing within a gel matrix provides evidence of new solid and semi-solid phase at low 
temperature J3]. The additional enthalpy observed here could indicate new solid phases 
for 80C-B+aerosil samples. These results are only observed upon initial crystallization of 
a freshly dispersed 80CB+aerosil sample. Bulk 80CB behavior is eventually approached 
upon repeated thermal cycling through the crystallization transition. Clearly, at these low 
silica densities, the silica gel is, at least locally, quite fragile and short range restructuring 
of the gel can strongly affect the liquid-crystal. 

The results of initial cycling through the crystallization transition for the ps = 0.220 



13 



sample shown in Fig. |H] do not reveal any significant changes. Aside from a very small, 
sharp, additional heat capacity feature at the iV-Sm^4 transition, the AC P curves are almost 
perfectly reproducible. The shift in Tin after the initial crystallization is small and upward 
by ~ 0.2 K, as expected by the expulsion of impurities. This indicates that for higher silica 
densities, the gel is robust and well behaved. Note that there is no appreciable change in 
the nematic range, as seen in both Fig. [7|and|Hl 



IV. DISCUSSIONS AND CONCLUSIONS 



Results have been presented from high-resolution ac-calorimetric experiments on 
80CB-t-aerosil dispersions with emphasis on the weakly first-order I-N phase transition. 
These results for 80C-B+aerosil dispersions have been compared with existing results for 
8C5+aerosil dispersions and reveal new aspects of the effect of quenched random 

disorder on liquid crystal phase transitions. In particular, these two LC+aerosil systems 
are very similar except for the relative elasticity of the LC material. The material 80CB is 
elastically stiffer, having an effective (single) nematic elastic constant larger by approx- 
imately 20%, than 8CB, and this is reflected by the higher transition temperatures for the 
nematic, smectic, and crystal phases. Thus, aspects that are LC material dependent and 
those that are general to quenched random disorder can be distinguished. 

From the very good overlap of the AC P wings away from the I-N two-phase coexistence 



region and from direct NMR studies on 8C-B+aerosil [Tjj a general feature of LC+aerosils 
is that the magnitude of the nematic order, S, is essentially the same as in the bulk LC. 
Thus, the main effect of QRD is on the director structure with the elasticity of the LC 
and the kinetics of the ordered phase growth as likely important factors. This is supported 
by the differences seen between 80CB+aerosil and 8CB+aerosil systems; in particular, 
the different ps dependence is likely connected to the difference in elasticity for the two 
liquid-crystals. 

The effect of quenched random disorder on first-order phase transitions is substantially 
different than that on continuous phase transitions. First-order transitions have an ad- 
ditional energy penalty for the formation of interfaces between coexisting phases, which 
complicates random- field type theoretical approaches. In the classical treatment, first de- 
veloped by Imry and Wortis 8], the QRD effect on first-order transitions is that a quenched 
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random-field creates domains having a randomly shifted transition temperature. This would 
have the effect of smearing the overall transition, a monotonic decrease in the transition la- 
tent heat and temperature with increasing QRD, and a low-temperature phase possessing 
only short-range order for arbitrarily weak QRD. These predictions are generally consistent 
with the behavior of nematic in aerogels 2^1 • F° r nematics in aerosils, the transition latent 
heat appears to decrease and the width of the coexistence region to increase monotonically 
with increasing disorder as well as the transition becoming apparently continuous for high 
disorder strength. These features are generally consistent with the classical picture. How- 
ever, the character of the transition and the non-monotonic transition temperature shifts do 
not appear to be consistent with this view. More strikingly, the I-N transition in aerosils 
in some range of ps, appears to proceed via two transitions. This could only occur in the 
classical view if a bimodal distribution of the random-field variance (h 2 ), connected to a 
bimodal distribution in p s , is present. This is not supported by SAXS studies, which re- 
vealed the fractal-like nature of the aerosil gel structure [lj , nor the behavior at the iV-Sm^4 
transition 3, However, nematics are very "soft" materials and the QRD imposed 

by aerosils appears to be much weaker than that of aerogels, thus elasticity of the LC (and 
possibly the gel) can play an important role. 

In LC+aerosil systems, (h 2 ) is thought to depend on the given concentration of silica and 
interaction with the LC, whereas the LC elasticity Kn is strongly temperature dependent 
(being proportional to S 2 ) near and below Tin- Since the aerosil gel is thixotropic and formed 
in the isotropic phase, any high energy strains or deformations that may exist are likely 
quickly annealed, the anisotropy of the gel should be fixed and essentially zero, especially 
for rigid gel structures like aerogels. However, the disordering nature of the gels may evolve 
with thermal history of the LC in aerosil gels. As seen by the result of cycl ing through 



crystallization presented here as well as the DNMR [17 1 and electro-optical jlG, 11] studies, 
the gel can be compliant with respect to distortions in the director structure for a range of 
silica densities. Note that the quantity Tin differs by at most 10% between 8CB and 80CB 
thus the energy scales are similar while the twist elastic constant differs by 36%. 

The features described above suggests a possible physical scenario for the origin of the 
double AC p feature at the I-N (or any "soft" first-order [31]) transition in LC+aerosil 
systems. As the nematic elastic constant strongly increases with decreasing temperature for 
T < Tin, a "skin" of low nematic order (due to the undulations of the aerosil strands) may 
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coat the silica strands. The thickness of this paranematic boundary layer would be strongly 
temperature dependent, shrinking with decreasing T below Tin- The presence of such layer 
a would serve to partially decouple the disordering (or field) effect of the silica gel from the 
void nematic (acting as a kind of "lubricant"). Once the layer thickness reaches its minimum 
value (roughly equivalent to a molecular length) the elasticity of the void nematic becomes 
strongly coupled to that of the aerosil gel. This would effectively increase, for a given ps, 
the disorder strength. 

A consequence of this speculation for a first-order transition induced change in disorder 
strength (through the onset of coupling between the director fluctuations and gel) would be 
the alteration of the gel dynamics (i.e., vibrational modes). This would be consistent with 
large changes in the relaxation times of aerosil gels observed by dynamic x-ray studies on 
8C5+aerosil near Tin [16]. In addition, this coupling should dampen director fluctuations 
and could account for the variation in the critical behavior seen at the iV-Smy4 transition 
for LC+aerosil samples. Another consequence of this view is that the director correlation 
length £ft (the relevant aspect of nematic order) would jump to a large isotropic value at 
the first transition and upon further cooling cross a second transition into a more strongly 
disordered state having a smaller correlation length. The isotropic nature of ^ and the 
final SRO state of the nematic with QRD are consistent with recent optical studies |ll( and 
the detailed evolution of £^ through the two transitions is the subject of current optical and 
calorimetric study j^]. 

The different silica density dependence of Tin and the temperature distance between 
the two C p peaks {ST2 P ) between 80C-B+aerosil and 8C-B+aerosil would also be consistent 
with the difference in the nematic elasticity of the two liquid-crystals. A stiffer silica gel 
(higher ps) would be required to influence a stiffer LC, thereby stretching the shift in Tin 
with respect to p$ as seen between the 80C-B+aerosil and 8CB+aerosil systems. The ps 
dependence in the high-density regime of 5T 2p for 80CS+aerosil (see Table HJ) compared to 
the nearly constant &Ti v ~ 0.1 K seen only in the low density regime of 8 CB+ aerosil 
would be compatible with a transition induced increase in the variance of the gel disorder 
strength. The liquid-crystal 8CB being much softer would only be able to stress a very 
weak (low p$) gel while the much stiffer 80CB would be able to distort a wider range of 
gels. Since the ps dependence of the critical behavior for the 7V-Sm^4 transition is quite 
similar between 8CB+aerosil and 80C-B+aerosil [4], these observations suggest that the 
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effects observed at the I-N transition are not directly connected to those at the iV-SmA 

Finally, the unexpected behavior of low silica density sample when initially cycled through 
the crystallization transition is not fully understood. The expulsion of silica impurities 
by the strongly first-order crystallization transition seems to lead to a more disordered 
system. One possibility to explain this phenomena is that the initial crystallization causes 
the expulsion of the silica particles locally and transforms the flexible fractal structure into a 
more rigid, foam-like, gel. The depression of T IN is then a consequence of the greater elastic 
distortions imposed by the new gel arrangement. Repeated cycling through crystallization 
would continue to expel impurities, eventually destroying the foam-like structure. Once 
the silica has been compacted sufficiently, percolation is no longer possible and free-floating 
silica particles would represent an annealed disorder. The gels formed in the high silica 
density samples are more robust and would not be expected to change significantly when 
the LC crystallize. This view can be directly tested with a detailed structural study by 
small-angle x-ray scattering where the thermal history is carefully controlled. The increase 
in the enthalpy of the iV-Sm^l transition after the initial crystallization to a value greater 
than the bulk LC value remains a puzzle. Note that the speculations presented here for the 
double I-N transition peaks and the unusual hysteresis behavior are intended to motivate 
future experimental and theoretical studies. 
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FIG. 1: Excess specific heat, AC P , as a function of temperature about Tjjv for bulk 80CB and 
80CI?+aerosil samples from pg = 0.036 to 0.647 grams of silica per cm 3 of liquid crystal. See 
figure inset for definition of symbols. The vertical dashed lines indicate the I-N transition region 
expanded in Fig. |2 
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FIG. 2: Expanded view of the excess specific heat about the I-N transition as a function of 
temperature. See figure insets for definition of symbols. The samples have been separated into 
two groups; the upper panel appears to indicate inhomogeneity induced broadening for samples 
with ps < 0.1; the lower panel depicts the evolution of AC P for ps > 0.1 and shows two distinct 
features, one sharp and one broad feature consistent with those seen in 8Ci?+aerosil systems 
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FIG. 3: The real (solid symbol, left axis) and imaginary (open symbol, right axis) I-N transition 
enthalpy are shown as a function of pg- The effective enthalpy 5Hj N is weakly dependent on silica 
content, indicating only minor changes to the latent heat conversion dynamics occur relative to 
the ac-frequencies employed in this work. The monotonic decrease in the imaginary component 
is evidence that the first-order character of the I-N transition decreases with increasing quenched 
disorder. Solid lines are guides to the eye while the horizontal dotted line represents the total 
pretransitional enthalpy 5Hjn = 5.13 J g _1 , which is independent of ps- See text for details. 
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FIG. 4: Dependence on ps of the I-N transition temperature, Tin, for 80C-B+aerosil (solid circles 
and left axis, T° N = 352.47 K) and 8CI?+aerosil (open circles and right axis, T° N = 313.99 K) 
samples. Data for 8C-B+aerosil samples taken from Ref. Note that both the left and right axes 
span 3 K in temperature. The solid lines are guides to the eye. 
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FIG. 5: Dependence on ps of the pseudo-transition iV-SmA temperature, T* scaled by the bulk 
value for 80C-B+aerosil (solid circles, T^ A = 339.52 K) and 8C-B+aerosil (open circles, T^ A = 
306.97 K) samples. Data for 8C-B+aerosil samples taken from Ref. if. The solid lines are guides 
to the eye. 
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FIG. 6: The nematic phase temperature range ATjv = Tin — T* scaled by the bulk value for 
80C£+aerosil (solid circles, AT N = 12.95 K) and 8C5+aerosil (open circles, AT N = 7.01 K) 
samples. Data for 8Ci?+aerosil taken from Ref. The solid lines are guides to the eye. 
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FIG. 7: Behavior of the I-N and iV-Sim4 excess specific heat of the ps = 0.051 sample as a 
function of temperature relative to Tjn before (solid circles, Tjn = 352.76 K ) and after (open 
circles, Tjn = 352.07 K ) sample crystallization. Both are heating scans made under identical ac- 
calorimetry conditions. Note the excess enthalpy for the AT-ShlA transition and the second feature 
near the I-N transition as well as a ~ 0.7 K shift downward of Tin observed after crystallization. 
See text for details. 
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FIG. 8: Behavior of the I-N and iV-Sim4 excess specific heat of the ps = 0.220 sample as a 
function of temperature about Tin before (solid circles, Tin = 351.02 K ) and after (open circles, 
Tin = 351.23 K ) sample crystallization. The after scan used the same ac-input power but a faster 
temperature scan rate than the scan before crystallization; both are heating scans. Note the nearly 
perfect reproducibility of AC P with a Tin shift upward of ~ 0.2 K observed after crystallization. 
See text for details. 
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TABLE I: Summary of the calorimetric results for 80C-B+aerosil samples. Shown are the sample 
density (ps in grams of aerosil per cm 3 of 80CB) as well as the I-N (Tin) and the TV-Sm^ 
(T*) phase transition temperatures, the nematic range (AT/v = Tin — T*), the width of the I-N 
coexistence region (STin), and the difference in temperature between the two AC P (IN) peaks 
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